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Abstract 

We investigate a mechanical system consisting of infinite number 
of harmonically coupled pendulums which can impact on two rigid 
rods. Because of gravitational force the system has two degenerate 
ground states. The related topological kink - likely the simplest one 
presented in literature so far - is a compacton, that is it has strictly 
finite extension. In the present paper we elucidate the relation of such 
system with sine-Gordon model. Also, solutions describing waves with 
large amplitude, and an asymptotic formula for the width of the kink 
are obtained. 
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1 Introduction 



Theoretical studies of classical mechanical systems with large number of de- 
grees of freedom are relevant to many branches of condensed matter physics 
(not to mention engineering sciences), where quantal aspects of microscopic 
world are practically not visible. Moreover, since the seminal work of Fermi, 
Pasta and Ulam [l] it is clear that such systems can also provide useful mod- 
els on which concepts of theoretical physics can be tested or illustrated. For 
example, studies of anharmonic lattices give insights into nonlinear transport 
phenomena, see e.g., (2], or mechanisms of thermal conductivity, see, e.g., jS]. 
Another example: sine-Gordon solitons, which appear in many places includ- 
ing quantum field theory (the duality with Thirring model), have beautiful 
realisation in a system composed of harmonically coupled pendulums, see, 
e.g., H. 

In recent papers jSt ifij we have investigated a mechanical system - de- 
scribed below in Section 2 - which has infinite number of degrees of freedom. 
It exhibits spontaneous symmetry breaking and related to it topological kink 
and antikink. Its main attractive feature is extraordinary simplicity of the 
kink. It has strictly finite extension equal to vr in a certain limit, and in 
that limit the kink is represented by the piece of sine function in the inter- 
val [— vr/2, 7r/2]. To the best of our knowledge, it is the simplest topological 
defect presented in literature so far. Furthermore, time evolution of the sys- 
tem during the symmetry breaking transition can be described analytically. 
For these reasons, the system is quite interesting as a testing ground for 
various ideas of theory of topological defects. Physics of such objects poses 
many very interesting albeit hard to solve dynamical questions, see [3 |H] for 
overviews. 

The purpose of the present paper is to supplement our previous investi- 
gations [Slini- In next Section we briefly recall the model and certain results 
in order to make the paper self-contained. We also introduce a folding trans- 
formation which maps trajectories of certain auxilliary, simpler model into 
trajectories of our system. Section 3 is devoted to large amplitude waves - 
the ones having small amplitudes were investigated in j^. The two types of 
waves involve very different kinds of motion. In Section 4 we calculate the 
width of the kink in a limit when it increases to inflnity. Connection of our 
system with the one considered in connection with the sine-Gordon soliton 
is elucidated in Section 5. We show that from mathematical viewpoint the 
unfolded version of our model can be regarded as non-analytically perturbed 
sine-Gordon system. In particular, the sine-Gordon soliton and breather pro- 
vide approximate solutions of our model. Nevertheless, our system and the 
sine-Gordon system are very different from physical viewpoint. Finally, in 



2 



Section 6 we point out several dynamical problems which can be posed within 
our model and which, in our opinion, deserve a careful investigations. 



2 The model and the folding transformation 

We consider one dimensional model with the Lagrangian 

L = \{dAf-\{d,ct>f-v{4>\ (1) 

where the field potential V{(t)) has the form 

1 COS0-1 for 101 <^o 
\ oo for 101 > 00, 

see Fig.l. Here 0o is a constant, tt > 0o > 0. Thus, all values of the real field 
are restricted to the interval [— 0o, 0o]- The potential has two symmetric 
minima at = ±0o, and it is not continuous at them. 





V(^) 

















Fig.l. The potential ^(0) 

This model describes an infinite set of harmonically coupled pendulums in a 
homogeneous gravitational field: ^s a dimensionless coordinate along a wire 
to which the pendulums are attached, r is a dimensionless time variable, and 
0(^, r) is equal to the angle between the upward vertical direction and the 
arm of the pendulum located at the point ^ and time r. The pendulums can 
move only in planes perpendicular to the wire. The restriction 

- 00 < </> < 00 (3) 
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is enforced by two rigid rods parallel to the wire. We assume that the pen- 
dulums bounce elastically from the rods. For a more detailed description of 
this system see |5j. 

Time evolution of r) is governed by the following equation 

(92-9|)0 = sin0 if 101 < 00, (4) 

and by the elastic bouncing condition 

5r0(e,r)->-9.0(e,r) if 0(e,r) = ±0o. (5) 

Hence, the time evolution of the velocities drcj) is in general not continuous. 
The vertical upward position, that is = 0, is unstable. The stable positions 
are given by = ±0o. Energy of the system is conserved. Typical motion of 
a pendulum consists of smooth 'flights' bewteen impacts with the rigid rods. 

In the case 0o ^ 1 we may replace sin0 by 0, and istead of Eq.(4) we 
then have 

(^2-91)0 = if 101 < 00. (6) 

Physically, the kink is obtained when the pendulums gradually turn from 
one rod to the other one along the wire. Its center is identified with the point 
^1 such that 0(^i) = - the corresponding pendulum points in the upward 
direction. Static kink is represented by a time-independent solution 0c(O of 
Eq. (4) or Eq.(6) which interpolates between — 0o and +0o, and such that 
0c and (9^0c are continuous when 0c — > ±0o. Such solutions have particularly 
simple form in the case of Eq. (6): for example, the kink located at = 
is given just by 

r -00 for e<-f, 
0c(O= 00 sine for -f<e<f, (7) 
I 00 for e>f- 

Anti-kink is represented by 0c(O = ~0c(O- Translations or 'Lorentz' boosts 
give, respectively, shifted or moving kink or anti-kink. As discussed in [H], 
the kink has strictly finite extension for any 0o < tt. Therefore, it is a simple 
example of so called compactons which were found for modified Korteweg- 
deVries models [HI [El E], and in a system of pendulums with non-trivial 
anharmonic coupling between them |12| . 

Condition (5) is quite troublesome because it introduces discontinuity in 
phase space trajectories of pendulums. For this reason, often it is convenient 
to 'unfold' the model. By this we mean passing to a new model with a 
new field 0(e, r) such that drcj) has continuous time evolution. can take 
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arbitrary real values. The relation between and has the following form. 
If 

^e[-(f)o + 2k(f)oAo + '2k(f)o\, (8) 
where k is an integer, then 

A(f: I ^(^' ^) ~ 2A;(/>o if k is even , . 

' ~ \ 2k(l)o-l(^,T) if A; is odd. ^ ^ 

The impacts on the rods occur when = 0o + 2Z0o with integer I. Relation 
(9) is depicted in Fig. 2. 




Fig. 2. The folding transformation 



Evolution of (j) is governed by the equation 

(9?-a|M=-^. (10) 

where V_{4>) is given by the following formula: if 4> belongs to the interval (8) , 
then in the case of Eq.(4) 

V{^) = cos(0 - 2k(f)o) - 1, 
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see Fig. 3, 



V 




Fig.3. The potential V{^) 



and 

H0) = -^(0-2A;0of (11) 

in the case of Eq.(6). 

Time evolution of (j) and of its first derivatives dr(p, d^cp is continuous. It 
follows from transformation (9) that condition (5) is automatically satisfied. 



3 Waves of impacting pendulums 

Let us first consider uniform motions of the pendulums, i.e., the ones such 
that (f) does not depend on ^. In this case Eq.(6) is reduced to 



if 

101 < 00- 
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The energy integral of motion corresponding to Eq.(12) has the form 



(t>'-r^ Co. (13) 

In the case the constant Cq is positive we may write Cq — v?, where > is 
the absolute value of the velocity of pendulums at the vertical upward posi- 
tion (0 = 0). This class of trajectories describes the pendulums impacting 
simultaneously on one barrier rod with the velocity {v? + 0o)^^^) elastically 
bouncing back from it and moving over the upward positon and further, un- 
til they impact on the other rod, from which they are reflected again - the 
motion is periodic in time. The pertinent solutions of Eq.(12) have the form 

0(t) = ±Msinh(r - tq). (14) 



At the time tq all pendulums are in the vertical position. They fall onto one 

ich that 

usvah{Tq - To) = 00. (15) 



of the rods at the time Tq such that 



Hence, 



T, - To = Arsinh(0o/^i) = In + ^1 + ^ j . (16) 

The period To of this motion is given by the formula 

ro = 4(T,-To). (17) 

The phase space picture of trajectories of this type is presented in Fig. 4 
by the curve A. Because the velocity changes its sign at the moment of the 
impact, the trajectory consist of two disconnected parts. 

The curves B and C in Fig. 4 represent two different sets of trajectories 
such that the pendulums do not pass through the vertical position (0 = 0): 
the pendulums bounce from one rod, raise a little bit and fall on the same 
rod. Also this motion is periodic in time. The corresponding solutions of 
Eq.(12) have the form 

0(t) = ±0^ cosh(r - To), (18) 



m 



where 0m is a constant such that < 0^ < 0o- In this case, co = —0^ 
in formula (13). At the time Tq pendulums reach their highest position, 
0(''"o) = =t0m- They fall on the rod after the time interval Tq — Tq, where 

cosh(T5 - To) = (19) 
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Hence, 



To 











I 0m \ 





(20) 



The period is equal to Tq 
were discussed in [S]. 



4(rg — To). Waves generated from these solutions 













By 


-<> o 













Fig. 4. Phase space trajectories of the uniformly impacting pendulums 
There are solutions of Eq.(12) for which Cq = 0: 

0±(r) = A0oexp[±(r-ro)], (21) 

where A = ±1. The solution 0+ describes pendulums which at the time tq 
start to move up from one of the rods with the initial velocity ^400 and reach 
the upward position at r = +oo. The other solution describes the reverse 
process. These motions are not periodic (Tq = oo), and they separate the 
previous three classes of trajectories. In Fig. 4 they are represented by the 
dashed lines. 

Finally, the thick dot in the center of Fig. 4 denotes the trivial, unstable 
trajectory = 0, and the thick dots at = ±0o the stable equilibrium 
positions. 

The spatially homogeneous solutions discussed above can easily be gen- 
eralized to solutions describing periodic in time and space waves just by 
performing a Lorentz-like boost. For example, in the familiar case of Klein- 
Gordon wave equation for a scalar field ip, = "^o^; instead of Eq.(12) 
we have (Pip/dr'^ = —m^il) with the solutions if) = A exp{±imoT) . Lorentz 
boost with velocity w, where \w\ < 1, amounts to the substitution r — > ^ = 
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7(t — w^), where 7 = (1 — w'^y^^'^. In this way we obtain the plane wave 
with the frequency on = ±7mo and the wave vector k = ±mo7w such that 
— k'^ = ml. 

Applying the boost to the solutions of the class S or C we obtain waves 
which coincide with the waves discussed in our previous work jB| after the 
substitution w = 1/v. Waves obtained from solution A are combined piece- 
wise from the functions 

0(r,O = ±Msinh(C-ro), (22) 

see Fig. 5. 




Fig. 5. The wave of impacting pendulums 

Here 

C = ^. (23) 

y/v^ — 1 

To is a constant, and v is the phase velocity of the wave. It is assumed that 
|f I > 1. The group velocity is equal to the boost velocity w = 1/v. The wave 
length and the frequency are equal to 

y 1 

Ao = 4:VV^ - l{Tg - To), ujo = —== — , 

— 1 -^0 
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where Tq is given by formula (17). Depending on the value of f , Aq can be 
any positive number, while < \uJo\ < 1/Tq. When f — > cxd we recover the 
uniformly bouncing pendulums: Aq = oo, coq = Tq^. 

The trick with boost can be applied if we have a solution which exists for 
all T - otherwise we would obtain 0(^, r) on a part of the ^ axis only. For 
this reason, solutions (21) have to be combined together in order to cover 
the whole time axis. There are two possibilities which differ by overall sign 
factor: 

^ / </'oexp(r - To) for r < tq, 
[ 0oexp(ro-r) for r > Tq, 

and — 0(t). Boosting these solutions consists of replacing r by C as above. 
We obtain a travelling wave in which the pendulum located at hits the 
rod at the time 

Tl = To H , 

V V 

while all other pendulums either are raising (those at ^ < ■Ci) or falling down 
(those at ^ > ^i). 



4 The width of the compacton when 0o ^r— 

In the particular case when (po = n the two rods merge into one put below the 
wire which supports the pendulums, and the angle <p can take arbitrary value 
in the interval [— 0o,0o]- As far as the static configurations are considered, 
such a system of pendulums coincides with the one corresponding to the 
sine-Gordon equation, where the rod is absent altogether. The point is that 
the least energy configurations of the systems coincide. Of course, for time- 
dependent configurations the two systems remain different because in our 
system the pendulums bounce from the rod, while in the sine-Gordon one 
they can move without encountering any obstacles. 

The soliton of sine-Gordon model has infinite extension because of its 
exponential tails. On the other hand, the compacton does not have tails for 
any 0o < tt. Therefore, its width ^o(0o) should be divergent when (j)^ ^ n 
from below. In [S^ we have obtained the following formula for the width 

eo(0o) = ^ /' d\ , \ (24) 

VZJo ^cos(A0o) - COS0O 

Substituting here 0o = vr — e, e > and changing the integration variable to 
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a = {tc — e)A/2 we obtain 



U<Po) = / da (25) 

^0 Jcos2(e/2) - sin^fT 



Second change of integration variable, namely 

sin 0" 



a ^ s 



cos(e/2) 

yields the elliptic integral of the first kind 



Uh) = t ds , ^ (26) 

;i-s2cos2(e/2))(l -s2) 







Thus, 



Because 



eo(vr-e) = 7r(cos^|). (27) 

1 16 
K{x) = - In- 



2 1 

for X —>■ 1—, see e.g. we finally obtain that 



eo(7r-e) = ln^ (28) 
sm 
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when e — ^ + . Thus, the width of the compacton diverges logarithmically 
when (f)Q —>■ TV—. 



5 Equivalence with the non-analytically per- 
turbed sine-Gordon system 

The potential ¥_{(/)) is periodic in (p, hence it can be written as Fourier series. 
Straightforward calculations show that potential (11) can be written in the 
form 

12 ^ (-1)" nTT 



m = -f 



■ cos 



(29) 



It is clear from Fig. 3 that the derivative V_ with respect to (j) does not exist 
at the points = fc0o, where k is an odd integer. At these points the Fourier 
series for dV_{<P) / dcf) is divergent. 
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The constant and the n = 1 terms in formula (29) give the potential 

ya)(0)^_^ + Meos(f0), (30) 
which is analytic in 0. Let us introduce the new field 

V'(V2T,v^e) = rf(e,T) + 7r. (31) 

Euler - Lagrange equation obtained from the potential K^^H^) t)e written 
in the form 

(a2-a|)z^(r,0 = -sin^(r,0, (32) 

which is identical with the sine-Gordon equation (on both sides we have 
replaced \/2t, \/2^ by r, ^) . 

The n > 1 terms in formula (29) give the non-analytic perturbation of 
the sine-Gordon potential (30), 

oo -1 

and of the sine-Gordon equation 

/P2 o2N / • / sin(n^) 

The sum on the r.h.s. of the last equation is finite for all values of ip, but 
it is not continuous at -0 = 2fc7r, where k is integer. This follows from the 
formula 

^ sin(n0) ^ , / • / 

n=2 ^ ^ 

and the fact that In function has a cut along the negative part of the real 
axis. 

The replacement of y(0) by y_^^\(f)) gives the sine-Gordon approximation 
to our system. The folding transformation applied to ^^''^^(0) yields the 
potential 

D TT^ 00 

where — 0o < < 00- It is presented in Fig. 6. The vertical lines at = ±0o 
in Fig. 6 correspond to the two rods. The main difference between this 
potential and the original one, see Fig. 1, is that now the first derivative of 
the potential vanishes at each of the two minima. This fact has profound 
influence on the shape of the soliton at large ^: the soliton acquires the 
exponential tails. 
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(1) 

V w 




Fig.6. The potential V^^^cj)). 



It is clear that with the help of the sine-Gordon approximation and for- 
mula (9) we can utilise the well-known exact solutions of sine-Gordon equa- 
tion in order to obtain approximate solutions of our system. For example, 
the sine-Gordon breather gives 




where < p < cxd is a parameter. The folding transformation transforms 
this solution into the approximate trajectory of our system of pendulums: at 
the time r = all pendulums are just bouncing from the rod at = —(f)o, 
next they move up. Some of them can cross the upward vertical position 
provided that p < I, but they do not reach the other rod. All pendulums 
stop simultaneously at the time 
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and they begin to slide back to the rod, from which they bounce again at the 
time 2Tjn-- The motion is periodic in time. We do not know whether there 
exists breather in our original system with evolution equation (4) or (6). 



6 Remarks 

1. Our model has several attractive features: simplicity of differential equa- 
tion (6); the simple analytical form (7) of the kink; neat physical realisation 
as the system of pendulums with their motions restricted by the two rigid 
rods. For these reasons we believe that the model can be useful for case 
studies of such interesting yet difficult to understand in detail phenomena as 
production of topological defects during symmetry breaking transitions (cer- 
tain preliminary results in this direction were presented in |n|), scattering 
of kink on anti-kink, or interaction of kink with a boundary in analogy to 
investigations reported in ^1^. 

2. Mechanical systems with impact have been investigated recently in con- 
nection with so called grazing bifurcation as a road to chaotic behaviour, 
see, e. g., [El CEl - the dynamics of such systems is far from trivial. Our 
system is of field theoretical type, that is it has infinite number of degrees 
of freedom, while most studies presented in literature so far are devoted to 
finite systems. It would be very interesting to find out whether the chaotic 
motions are present in our system and, if present, what is their role in dy- 
namics of compactons. 
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